Coupled-cluster theory for three-body Hamiltonians 
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We derive coupled-cluster equations for three-body Hamiltonians. The equations for the one- and 
two-body cluster amplitudes are presented in a factorized form that leads to an efficient numerical 
implementation. We employ low-momentum two- and three-nucleon interactions and calculate the 
binding energy of *He. The results show that the main contribution of the three-nucleon interaction 
stems from its density-dependent zero-, one-, and two-body terms that result from the normal 
ordering of the Hamiltonian in coupled-cluster theory. The residual three-body terms that remain 
after normal ordering can be neglected. 
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I. INTRODUCTION 

One of the central challenges in nuclear theory is to 
understand and predict the structure of nucleonic matter 
based on microscopic nucleon-nucleon (NN) and many- 
nucleon interactions. In recent years, there has been sig- 
nificant progress in exact calculations of ground and ex- 
cited states of light nuclei based on various high-precision 
interactions fitted to NN data [H, 0, H, 0, H, 11, H] . These 
results clearly show that three-nucleon forces (3NFs) con- 
tribute significantly: Without 3NFs, the binding energies 
depend strongly on the NN potential used, which can 
be traced to scheme and model dependences in any the- 
ory restricted to NN interactions. The existence of 3NFs 
is not surprising, since nucleons are not point particles. 
There are always virtual excitations (high-momentum 
nucleons and A-isobars) or internal degrees of freedom 
(quarks and gluons) that have been "integrated out" . 
This directly leads to three- and many-nucleon inter- 
actions. The modern understanding is that nuclear in- 
teractions are effective interactions and depend on the 
resolution scale given by a cutoff A of the effective the- 
ory (see for example Q). Exact calculations are cutoff- 
independent up to the effects of omitted higher-order in- 
teractions, and therefore varying the resolution scale is a 
powerful tool to analyze the predictive power of theoret- 
ical calculations. 

The study of SNFs in systems beyond the lightest nu- 
clei is an important goal. This requires a flexible tech- 
nique to solve the many-body problem including NN and 
3N interactions. Coupled-cluster theory is a promising 
tool for this endeavor. This method originated in nuclear 
physics [1, [13] and is today mostly propell ed through its 
importance in quantum chemistry llll . 11211 . For reviews, 
we refer the reader to Refs. After the 

seminal work by the Bochum group (l3| , Heisenberg and 
Mihaila employed coupled-cluster theory for structure 
calculations of ^^O based on realistic NN potentials p^ . 



For their calculation of the charge form factor of ^^O, 
they also included selected contributions from 3NFs that 
could be cast into the form of density-dependent NN in- 
teractions [T9I. Another recent approach employed ab- 
initio coupled-cluster theory for structure calculations in 
closed-shell nuclei ^He and i^Ojp,|2l|,l2i|, in open-sheU 
nuclei as the neighbors of ^^O [23|, and weakly bound and 
unbound helium isotopes [13]. These calculations were 
limited to NN interactions. It is the purpose of this paper 
to develop coupled-cluster theory for three-body Hamil- 
tonians. This extension of the coupled-cluster method 
might also find applications in condensed-matter theory 
and quantum chemistry [2^, [2I, |27j . 

A current frontier in nuclear structure theory is to 
determine consistent 3NFs corresponding to the differ- 
ent NN interactions, and with predictive power when 
extrapolated to the extremes of isospin and to moder- 
ate densities. Several theoretical approaches are cur- 
rently being used. The Tucson-Melbourne 3NF al- 
ready employed symmetries of QCD in its construc- 
tion [2^ [2^. Existing phenomenological 3NFs include 
the Fujita-Miyazawa force based on 27r exchange with 
an intermediate A-isobar [30| and also various shorter- 
ranged 3NFs m, m, I34j . This approach has led to 
a very successful description of light nuclei (for a review 
see Ref. 0). 

Onother approach is to systematically construct NN 
and higher-order interactions within the framework of 
chiral effective field theory (EFT) [H, H 113, HI . Chiral 
interactions are expanded in powers of a typical momen- 
tum of nucleons in nuclei and in powers of the pion mass, 
both generically called Q, over the EFT breakdown scale 
A^ ~ 1 GeV (« 5fm~^). The EFT power counting natu- 
rally explains the hierarchy of NN, 3N, and higher-body 
interactions [SO^, which only enter in subleading orders 
and make calculations for complex nuclei based on NN 
and 3N interactions meaningful. At this point, the lead- 
ing 3NF has been implemented [iO, Hi], \^ . 
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Nuclear interactions require regularization and renor- 
malizaton to be meaningful, and with EFT, usually a 
momentum cutoff scheme is used [s^, [H, Es, The 
description of the NN data is not completely independent 
of the cutoff, but the cutoff variation decreases with in- 
creasing order, because the cutoff dependences can be ab- 
sorbed by higher-order contact interactions. In this way, 
chiral EFT implements the renormalization group (RG) 
running of the interactions up to higher-order terms. 
This is reflected in the leading chiral 3NF, which includes 
a 3N contact term to be fitted to 3N data. In principle, 
nuclear structure calculations based directly on chiral in- 
teractions are feasible QjISIIIjIISI- However, the typical 
cutoffs (A « 2.5-3fm"^) are somewhat too large to use 
chiral interactions without resummations or prediagonal- 
izations in many-body calculations. 



Renormalization group methods can be used to evolve 
nuclear interactions to lower momenta, which leads to 
improved convergence in few- and many-body calcula- 
tions [i^ [50, The resulting low- momentum 
interactions, known generically as "Mowfe", have vari- 
able momentum cutoffs and are approximately indepen- 
dent of the starting NN interaction for cutoffs A < 
2fm"^ [MSi- The RG evolution preserves the long- 
range parts, and starting from chiral EFT interactions, 
generates all higher-order contact operators needed to re- 
produce low-energy NN observables. Moreover, with in- 
creasing orders in EFT, the resulting Mowfc interactions 
are very similar to low-momentum interactions obtained 
from conventional potentials [i^ [s^l • Since chiral EFT 
represents the most general low-momentum expansion of 
nuclear forces, the above observations motivate combin- 
ing low-momentum interactions with 3NFs from chiral 
EFT [5OI . We will follow this approach and employ Mow fc 
with the corresponding low-momentum 3NF adjusted to 
the binding energies of and *He [50|. In this way, we 
can expect to define approximately consistent NN and 
3N interactions. In this first study, we focus on the de- 
velopment of the coupled-cluster method to include 3NFs 
and will only present results for one cutoff. A study of 
the cutoff variation and associated uncertainties due to 
higher-order many-body interactions will be left to future 
work. 



This paper is organized as follows. In Section [IT] we 
present coupled-cluster theory and its extension to three- 
body Hamiltonians. In Section [1111 we solve the resulting 
coupled-cluster amplitude equations at the singles and 
doubles (CCSD) level for ^He based on low-momentum 
NN and 3N interactions. Our findings are especially 
promising as they show that the binding energy of "^He 
can be calculated based on those parts of the 3NF that 
can be viewed as density-dependent zero-, one-, and two- 
body forces. We summarize our results in Section HVl 



II. COUPLED-CLUSTER EQUATIONS FOR 
THREE-NUCLEON FORCES 

This is the main technical section of this paper. In the 
first subsection, we briefly recapitulate coupled-cluster 
theory and the reformulation of the three-body Hamil- 
tonian in normal-ordered form. The second subsection 
deals with the diagrammatic derivation and the factor- 
ization of the coupled-cluster equations due to the resid- 
ual three-body force that remains after normal ordering 
of the Hamiltonian. 



A. Coupled-cluster theory 

We consider a pure three-body Hamiltonian. Coupled- 
cluster theory for one- and two-body Hamiltonians is a 
mature field, and we refer the reader to the reviews [H, 
El, [H, [M [l3 ■ The three-body Hamiltonian is written 
as 

= ^ X! {P<ir\\stu)alalalauatas . (I) 

pqrstu 

Here, (pqr\\stu) denotes the antisymmetrized three-body 
matrix elements, while a| and cip create and annihilate a 
fermion in the single-particle orbital p, respectively. 

In coupled-cluster theory, the Fermi vacuum is a single- 
particle product state = HiLi^llO), where the A 
lowest-energy orbitals are occupied. In a first step, we 
cast the three-body interaction into a normal-ordered 
form with respect to this vacuum. In what follows it 
is assumed that the indices i, j, fc, Z, m label the occupied 
orbitals of \4>) while a,b,c,d,e refer to the unoccupied 
orbitals of \(f>). The former indices run over the number 
Uo = A of occupied orbitals, while the latter run over 
the remaining number rt„ of unoccupied orbitals. Typi- 
cally, one has 7i„ 3> rio- Indices referring to all orbitals 
are denoted as p, q, r, s, t, u; see for example Eq. ([T]). The 
normal-ordered Hamiltonian is thus 

ijk ijpq 

+ ^ X! {■ipq\\-irs){aldldsdr} + , 

ipqrs 

where /i3 denotes the residual three-body Hamiltonian 
= ^ X! {P1^\\^^^){°'l°'\°'l°'u^tds} ■ (3) 

pqrstu 

Here, we used the {. . .} to denote normal ordering. Writ- 
ing the three-body Hamiltonian in normal-ordered form, 
it is clear that the Hamiltonian separates into a zero-, 
one-, two-, and a three-body term. The first three sums 
in Eq. ^ are the vacuum expectation value, and the 
"density-dependent" one- and two-body terms, respec- 
tively. Their treatment is standard in coupled-cluster 



theory, as they simply modify the normal-ordered two- 
body Hamiltonian from the NN interaction. In this sec- 
tion, we focus on the residual three-body operator of 
Eq. ©. 

In coupled-cluster theory, the correlated state {ip) is 
given by a correlation operator exp (T) that acts onto 



a single-particle product state 
A-body system by 



nf=ia!lO) of the 



(4) 



The cluster operator, 



T = Ti+T2 + 



Ta, 



consists of a one-body cluster operator 



(5) 



(6) 



FIG. 1: (Color online) Energy contributions of the residual 
three-body Hamiltonian, Eq. ([3)l, in the CCSD approxima- 
tion. 



into the first of the CCSD equations, Eq. ([S]), to deter- 
mine the ground-state energy. Note that the similarity- 
transformed Hamiltonian of Eq. pT|) is not Hermitian, 
and CCSD is not a variational approach. However, the 
similarity-transformed Hamiltonian can be evaluated ex- 
actly for any truncation of the cluster operator. In what 
follows, we will compute the corrections to the energy, 
Eq. (O, and to the CCSD Eqs. dH) and that arise 
due to the residual three-body Hamiltonian. 



a two-body cluster operator 

^2 = 7E^^''«1«I«^- 



(7) 



ijab 



and so forth. Note that the two-body cluster ampH- 

_^6a ^ab 



tudes tff 



are fully antisymmetric. 



Clearly, Ti and T2 induce \p — Ih and 2p — 2h excita- 
tions, respectively. In what follows, we will limit the 
expansion, Eq. ([S]), of the cluster operator to the one- 
body cluster, Eq. ([6]), and the two-body cluster, Eq. ([7]), 
respectively. This approximation is referred to as CCSD 
("coupled-cluster theory with single and double excita- 
tions"). CCSD is a powerful approximation and a com- 
promise between accuracy on the one hand and compu- 
tational effort on the other hand. One inserts the ansatz 
Eq. ^ into the Schrodinger equation, multiplies with 
exp (— T) from the left, and obtains the following set of 
equations 



E = 
= 
= 



b\H\4>) 



Here IC^'.V.f^") 

citation of the product state \q 



and 



(8) 
(9) 
(10) 

is a — nh ex- 



H = cxp(-T)i/cxp (T) 



(11) 



is the similarity-transformed Hamiltonian. This Hamil- 
tonian is a sum of one-, two-, and three-body Hamilto- 
nians, i.e. H = Hi + H2 + Hs- The treatment of this 
Hamiltonian within coupled-cluster theory is well known, 
except for the residual three-body term, Eq ([3]), of the 
normal-ordered three-body Hamiltonian H3 . 

The CCSD Eqs. dH) and ^ determine the one- 
particle and two-particle cluster amplitudes tf and t^j", 
respectively. These amplitudes can then be inserted 



B. Derivation of coupled-cluster equations 

In this subsection, we derive the contribution of the 
residual three-body term ^13 from Eq. ^ to the energy, 
Eq. dH]), and the cluster amplitudes, Eqs. ([9]) and (fTO| . 
The final results are given in Eq. (jl2p for the energy and 
in Eqs. (fTS]) and for the cluster amplitudes, respec- 
tively. The contributions to the energy and the coupled- 
cluster amplitudes can be derived, for instance, by follow- 
ing the approach of Ref. [lB| • The resulting expressions 
contain 2, 15 and 51 terms for the energy and the clus- 
ter amplitudes. Many of the individual terms consist of 
sub-terms of similar structure. The theoretical deriva- 
tion presented in this subsection exploits this structure 
and leads to an efficient numerical implementation. 

It is most convenient to evaluate the matrix ele- 
ments in Eq. ([8]) of the similarity-transformed Hamilto- 
nian. E g. (TT|) . in a diagrammatic form (see, for exam- 
ple, [lj,l53l). Dia grams are a useful book keeping device 
to keep track of the (considerable) number of possible 
Wick contractions. We refer the reader to the literature 
for a more detailed description. 

In a first step, we determine the correction to the 
CCSD energy, Eq. ([8]), that is due to the residual three- 
body Hamiltonian The matrix element (0|ft,3|0) is 
a sum of all topologically different diagrams, where the 
Hamiltonian, Eq. ([3]), is fully contracted by the cluster 
operators. The two diagrams that enter this expression 
are presented in Fig. [TJ In these diagrams, the thick 
horizontal bar represents the residual three-body Hamil- 
tonian, Eq. ([3]); the thin horizontal bars denote the one- 
body and two-body cluster operators, respectively. Par- 
ticle and hole lines are also shown. If uncontracted, the 
former have an arrow that points upward while the latter 
have an arrow that points downward. The corresponding 
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+ 



07i\/ + n f^^A/ + (^VA/ 



FIG. 2: (Color online) Contributions of the three-body Hamil- 
tonian, Eq. ([3]), to the Ti cluster equation in the CCSD ap- 
proximation. 



algebraic expression for the energy correction is 

63 = {<P\{Hf^f^)^\<f>) + {^\{^Hfl)^\4,), (12) 



- J2 {klm\\cde)tltt 



klmcde 



i ^ {klm\\cde)tltft'^. 



klmcde 

This is the energy correction due to the residual three- 
body Hamiltonian. Note that the computational effort of 
the energy, Eq. (fT2|l . scales as 0{nl^nl). Note also that 
the summation over J2kc common to both diagrams 
and might therefore be factored out. This is the basic 
idea behind the factorization [s^l, and will be presented 
in detail below. 

Let us consider the matrix element {(j)'^\H\(j)) appearing 
in the CCSD Eq. ([9]). It is given by the sum over all 
topologically different diagrams where one particle line 
and one hole line are not contracted. Figure [2] shows the 
15 diagrams that enter this matrix element. 

The translation of these diagrams into algebraic ex- 
pressions is straightforward, but there are two reasons 
to not present them here. First, the naive numerical 
implementation of the resulting expressions would be in- 
efficient. To see this, consider, for instance, the third 
diagram in Fig. O The corresponding algebraic expres- 
sion is 



1 



{klm\\cdi)tlt 



da 
Im ' 



(13) 



ckdlr 



and its naive evaluation costs 0{n^n'^) operations. How- 
ever, performing the summations involving the Ti cluster 
operator first yields the intermediate 



jhn 



]-Y{klm\\cdi)tl , 



(14) 



which only costs 0{n^n'^) operations and requires 
0{nun^) in memory. The subsequent contraction of this 



intermediate with the remaining T2 cluster operator only 
costs ©(n^n^) operations. Clearly, the memory cost of 
the intermediate is overcompensated by the reduction of 
computational cycles. Note also that the intermediate 
Eq. (|14p enters in the evaluation of the seventh diagram 
depicted in Fig. [21 The second reason is that the com- 
plexity of the involved diagrams increases rapidly. The 
number of diagrams increases from 2 to 15 to 51, when 
going from Eq. ^ to Eq. ([9]) to Eq. (flO)) . respectively. 
The construction of each individual diagram and its in- 
spection regarding the construction and use of intermedi- 
ates then becomes cumbersome, and a more systematic 
approach is called for. Similar comments apply when 
improving the coupled-cluster wave function through the 
inclusion of three-body or four-body cluster amplitudes. 
One therefore considers a factorization of the coupled- 
cluster equations . This approach yields a very com- 
pact form of the coupled-cluster equations and is particu- 
larly useful for the numerical implementation [54, 55, 5^]. 
So far, factorized coupled-cluster equations have been de- 
rived in a two-step procedure. The first step consists of 
constructing all topologically different coupled-cluster di- 
agrams. In a second step, these diagrams are analyzed 
and repeatedly decomposed into simpler intermediates 
that undergo single contractions. Here, we proceed dif- 
ferently and present a direct diagrammatic derivation of 
the factorized coupled-cluster equations. Our derivation 
avoids the explicit construction of all individual coupled- 
cluster diagrams. 

We have to decide in which order multiple contractions 
of the Hamiltonian with the cluster operators should be 
performed. Let Tk{p"h™) denote the contraction of the 
Hamiltonian with the fc-body cluster operator Tk via 
n particle lines and m hole lines. The contraction of 
Tk{p"'h™) with an object of i particle lines and j hole 
lines costs cn\jil, computational operations and results 
in an object of size mn\ji^^. Here, c and m denote the 
computational cost and memory multiplier, respectively, 
and one finds 



Based on this analysis, we find that the cost of two sub- 
sequent contractions labeled Ta and Tb, respectively, is 
proportional to c(Ta) + •m{TA)c{TB) when contraction 
Ta is first, and proportional to c(Tb) + m{TB)c{TA) 
when contraction Tb is first. Table [1] shows (from left 
to right) the optimal order in which subsequent con- 
tractions should be performed, under the condition that 
1 ^ rio <C n„ and < n„. We also listed the relative 
computational cost c, and the memory multiplier m. For 
the first 7 entries, the order is easily understood. These 
contractions do not yield an increase of the size of the 
contracted object (since m ^ 1), and the order is there- 
fore determined by the computational cost. The remain- 
ing four contractions increase the size of the contracted 
object (since m > 1), and it is usually most efficient to 
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TABLE I: Memory multiplication factor m and computational cost factor c for forming intermediates through contraction with 
the fc-body cluster operator Tk{p"h^^) via n particle lines and m hole lines. 







Tiiph) 




Tiip) 


T2{p') 


T2{ph^) 


T2{ph) 


T2(p) 


T2(h^) 


Ti(h) 


T2{h) 


m 








Uo/Uu 






1 


nl 




Tin /no 




c 


1 


1 


rio 


rio 


nl 








nl 







p4 



□51 



+ 



+ 



OP 



:7V 



:7V 



VDT 



:ir 



IT 



1 EZb; 

2 iT^ 



IT 



FIG. 3: (Color online) Matrix elements of the residual 
three-body part of the similarity-transformed Hamiltonian, 
Eq. ([TT|l. that enter the amplitude Eq. ([9}. 



line. These diagrams should be constructed from simpler 
diagrams, adding one contraction at each step. For the 
residual three-body Hamiltonian, the root is clearly given 
by the diagram 1 in Fig. [3] This diagram is the only one 
that has three incoming particle lines and three outgoing 
hole lines. Diagrams 2 and 3 have one outgoing particle 
line and one incoming hole line, respectively, and also 
have a total of six outgoing and incoming lines. They 
are sums of two diagrams. The first is the residual three- 
body Hamiltonian with this appropriate particle and hole 
lines, while the second diagram is a contraction of dia- 
gram 1 with a Ti cluster operator. 



In what follows, we adopt the following convention. 
For a diagram that is the sum of diagrams, we label the 
first, second, third, etc. term of the sum by a, b, c, re- 
spectively. For example, diagram 2 is the sum of diagram 
2a and diagram 2b. These labels are not printed in Fig. [31 
Diagram 4 is the sum of 4 diagrams, namely the corre- 
sponding residual three-body Hamiltonian and three con- 
tractions of previously generated diagrams. Note that 
out of the two diagrams of diagram 2, only diagram 2a 
enters. This is due to the order specified in Table HI Di- 
agram 2b is a contraction of diagram 1 with Ti{p), and 
this contraction cannot be followed by the contraction 
Ti{p). It is now clear how to proceed. Diagrams 5 to 8 
have a total of four incoming and outgoing lines, while di- 
agrams 9 to 12 have a total of two incoming and outgoing 
lines, respectively. Note the peculiar factor 1/2 in front 
of diagram 9b. This diagram is a Ti{ph) contraction of 
diagram 5b which itself is also a Ti{ph) contraction. The 
factor 1/2 will be needed for the translation into alge- 
braic expressions. Similar comments apply to diagrams 
10c, 11c, and 12d. Note that diagram 12 consists of all 
diagrams with one incoming hole line and one outgoing 
particle line. The recursive expansion of the correspond- 
ing right-hand side yields indeed all diagrams depicted in 
Fig.m and therefore factors the CCSD Eq. ^. 



perform the computationally more expensive contraction 
before performing the second contraction on an object 
with increased size. 

Let us now turn to the diagrammatic factorization of 
the coupled-cluster equations. Considering Eq. ^ and 
Fig. [21 we have to construct all coupled-cluster diagrams 
with one incoming hole line and one outgoing particle 



Let us translate Fig. [2l into algebraic expressions. We 
use the convention that the intermediate Ipl'pl'"' has in- 
coming lines pi,p2, ■ ■ ■ and outgoing lines gi, g2, ■ • re- 
spectively. In what follows, the intermediate corre- 
sponds to diagram i/ of Fig. [3l We restrict ourselves to 
those intermediates that are needed for the construction 
of diagram 12, and denote the diagrams 1, 2a, 3a, and 4a 
directly in terms of the corresponding three-body matrix 



6 



TV" 



~2T" 

-TV 



t t 



TV 



s + 



TV" 



IT 



+ 



"IT 



t V 



+ 



7V 



+ j 



+ T 



Ho 



or 



♦ 



+ 



y W 

or 



FIG. 4: (Color online) Additional diagrams needed for the contributions of the residual three-body Hamiltonian, Eq. ([3|, to 
the CCSD Eq. ([TUJ. 



elements. The result is 



and 



d 



km 



e 



/(lOac): 
/(ll)r 



^J2{kla\\cde)tl'} + Y,I(6htX,, 

ckdl kc 

-Y.{klm\\cd^tfi+2Y^I{btTtl 



ckdl 



kce 



kc 



In terms of these intermediates, the correction to the 
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right hand side of the CCSD Eq. ^ reads 

^Y.{kla\\cdt)ti:l + J2l{6h)'!^Jl!l (15) 

ckdl eke 



and 



ekni 



kc 



ek 



E^(10ac):i^-E^(ll)rC 



These terms have to be added to the weU-known CCSD 
equation for the Ti-cluster amphtudes based on two-body 
Hamiltonians. 

For the factorization of the CCSD Eq. (|10p , one has to 
construct all diagrams with two incoming hole lines and 
two outgoing particle lines. Their number is 51, and we 
directly derive them in a factorized form. Figure |4] shows 
the factorized diagrams that are needed for the recursive 
construction of the CCSD Eq. (fTO)l . The corresponding 
algebraic expressions for the intermediates are 



/(3) 



klm 
cdi 



12 
1 

12 



{klm\\cd^)+I{3bt!iT, 



e 

/(7c)^r = 2E^(5)L"^^ 



\k7n 
)ei 1 



Ide 



Imd 

mt^ = 4E^(5)L"c., 

me 
e 

mabcett = /(8ac)^7 + /(86)^,^+/(8e)l^ 
I{lQabc)l = /(lOac): - 2E ^(5)L"C™ , 



kmc 



l-{klmm) + ^J2^klm\\cde)tl 

ce 

J2(l{klrn\\cdj)+I{3b)fjAt^,, 



I{^^)fe = \ Y.{alb\\cde)tt 

Id 

^(16)^'" = \T.'^klm\\cdj)tt+\Y.I{ht:^t% 

cdl ce 

+ E (^(7a&)cf + \wc)1t) n 

+ \ Y.^klm\\idj)tf , 

Id 

im-d- = i(fc/&iizdj) + iE^^^^ii^^^)*?i 

ce 

+ eQ(^^^im-?')+^(4c)c*)^^ 

/(19ac)^^ = IE^^^^II^^^)^"/ + E<"^^II^^-?')*'' 

/fie Id 

I{m)f^ = 2E^(15)cei|, 

e 

old ee 

+ 2E/(7rcf4^ + ^E(^^^iMj>' 

mc Id 

+ E(n8ac)^)+/(8e)^))i^ 

c 

+ ^E^(9)'^^i'-3E^(14)Sf* 

C 

1(20^^" = -2E^(i6)^'"t: 



db 
Im 7 



Idm 



The residual three-body Hamiltonian, Eq. thus 
leads to the following correction, to be added to the right 
hand side of the CCSD Eq. 
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Y,{alb\mtf + ip(zj) Y,{alb\\cdj)tfr + P{ab)Pitj) J2 K(8a6ce)^f + ^/(8d)^n tl"^ (16) 

Id led kc ^ ' 

e m ce 

P(a6)P(zj) E ^(18)^)4" + Pi^j) E (^(19«c)c' + ^^(19&)?}) 

feid c ^ 

P(a6)P(zj) E (^(20a - g, i)^ + ^1(20/^)^'') tl . 



km 
\ab \ ic 



Here the permutation P{ab) implies P{ab)Iab = lab^ha- 
Again, these have to be added to the right hand side 
of the CCSD Eq. Equation ^ for the energy 

correction, and the expressions of Eqs. p5)) and 
for the cluster amplitudes are the main technical re- 
sults of this paper. Upon expansion of diagram 22 in 
Fig. m one gets indeed all 51 diagrams that enter the 
coupled-cluster Eq. (jlO[) . The numerical implementa- 
tion of the terms in Eqs. (fT5| and p6|) is straightfor- 
ward. It is interesting to analyze the resulting compu- 
tational costs. The most expensive intermediates /(IB) 
and /(19ac) cost nf^n^ and n^nl computational cycles, 
respectively. The most memory-expensive object is the 
interaction {alb\\cde) which enters the construction of in- 
termediates /(15) and /(19ac), and requires the storage 
of n^Tio real numbers. In our largest calculations for ^He, 
we have n„ = 4 and n,, = 220. 
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III. APPLICATION TO ''He 

We present the first ab-initio coupled-cluster calcu- 
lations includin g 3N Fs. Our results are based on low- 
momentum NN |49j and 3N [s^] interactions, 

i/ = r + Mowfe(A) + y3N(A). (17) 

In this exploratory study we use a sharp cutoff A = 
1.9 fm^^, and Viowfc is derived from the Argonne vis po- 
tential |57j- The corresponding 3N interaction is based 
on the leading chiral 3NF and has been fitted to the •^H 
and ^He binding energies in Ref. [soj. This 3NF consists 
of a long-range 27r-exchange part, determined by the low- 
energy coefficients ci, C3, and C4, an intermediate-range 
iTT-exchange ( "D-term" ) and a short-range contact inter- 
action ( "i?-term" ) 1140', 42] . The 3NF operators are multi- 
plied by regulating functions of the incoming and outgo- 
ing Jacobi momenta, /r(p, q) = exp[-((p2 + 3^)/A2) ] , 
with the same cutoff value A as in T^ow k ■ For additional 
details on the 3NF, we refer the reader to Refs. [50ll5l|. 
These low-momentum interactions can be directly em- 
ployed within coupled-cluster theory. 

In this application of coupled-cluster theory we re- 
strict ourselves to a proof-of-principle calculation. For 



FIG. 5: (Color online) CCSD results for the binding energy 
of *He as a function of the oscillator spacing and for model 
spaces consisting of A = 3 to A = 6 oscillator shells. The 
CCSD calculations are based on low-momentum NN and 3N 
interactions, where the full and dashed lines respectively de- 
note the energy obtained with and without 3NFs. 



simplicity, we have therefore considered only the 3NF 
channel with total isospin T = 1/2, total angular mo- 
mentum J = 1/2, and positive parity. This partial 
wave is the dominant contribution to the binding ener- 
gies of light nuclei: For *He, the corresponding Faddeev- 
Yakubovsky resuh is £; = -28.20(5) MeV, which differs 
only by 100 keV from the exact energy including all par- 
tial waves, E = -28.30(5) MeV [13]. 

The coupled-cluster calculations are performed in a 
harmonic oscillator (HO) basis, with basis parameters 
given by the oscillator spacing hiu and the number N 
of oscillator shells. We will present our results as a 
function of these parameters. The matrix elements of 
the 3NF are calculated first in relative HO states given 
by the expansion in 3N partial waves. For the trans- 
formation from relative HO states to the single-particle 
("m-scheme") basis, we essentially follow the Appendix 
of Ref. fi\. However, we can also start from a 3NF in 
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FIG. 6: (Color online) Data points; CCSD results (taken at 
the fiw minima) for the binding energy of ''He with 3NFs as 
a function of the number of oscillator shells. Dashed lines: 
Exponential fit to the data and asymptote of the fit. Full 
line: Exact result. 



the non-antisymmetrized HO basis, which simplifies the 
transformations given in Eqs. (B.9), (B.ll), and (B.12) 
of Ref. and explicitly antisymmetrizc the m-scheme 
matrix elements in the last step of our transformation. 
This approach does not require coefficients of fractional 
parentage. We verified that our matrix elements agree 
with those obtained from a transformation based on an- 
tisymmetrized relative HO matrix elements. We further 
have checked that the transformation preserves the unit 
matrix and that it yields identical matrix elements for 
(Ti • (T2 and (J2 ■ T3 in antisymmctrized states. 

In Fig. [5l our CCSD results for the binding energy of 
"^He are shown with and without 3NFs as a function of 
the oscillator spacing and with increasing model space 
size. The 3NF contribution is repulsive, in agreement 
with the Faddeev-Yakubovsky calculation [HO], and cor- 
respondingly, the minima in Huj are shifted to smaller os- 
cillator spacings. We observe a slow convergence at the 
last few 100 keV level, which is due to the sharp cutoff in 
Mow k ■ This might be improved b y u sing low- momentum 
interactions with smooth cutoffs [58|. Using the minima 
of the CCSD results with 3NFs, we make an exponential 
fit of the form E{N) = Eao + aexp {—bN) to the data 
points. The result is shown in Fig. [6l The extrapolated 
infinite model space value is Eoo = — 28.09 MeV, which 
is very close to the exact result E = —28.20(5) MeV. 

It is interesting to analyze the different contributions 
AE to the binding energy E. The individual contribu- 
tions are given in Fig. [7] for a model space of TV = 4 
oscillator shells and hoj = 20 MeV. The main con- 
tribution stems from the low-momentum NN interac- 
tion. The contributions from SNFs account only for 
about 10% of the total binding energy. This is con- 
sistent with the chiral EFT power-counting estimate 
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FIG. 7: (Color online) Relative contributions \AE / E\ to the 
binding energy of *He at the CCSD level. The different points 
denote the contributions from (1) low-momentum NN inter- 
actions, (2) the vacuum expectation value of the 3NF, (3) the 
normal-ordered one-body Hamiltonian due to the 3NF, (4) 
the normal-ordered two-body Hamiltonian due to the 3NF, 
and (5) the residual 3NFs. The dotted line estimates the 
corrections due to omitted three-particle-three-hole clusters. 



(Msn) {Q/Ax)L^iov,k) « 0.1 (Mowfe) m (see also 
Table 1 in Ref. [52|). The second, third, and fourth 
largest contribution are due to the first, second, and 
third term on the right-hand side of Eq. ([2]). These 
are the density-dependent zero-, one-, and two-body 
terms, which resulted from the normal ordering of the 
three-body Hamiltonian in coupled-cluster theory. The 
contributions from the residual three-body Hamiltonian, 
Eq. ([3]), are very small and are represented by the last 
point in Fig. [T] Recall that the residual 3NF contributes 
to the energy directly through Eq. p2|) and indirectly 
through a modification of the cluster amplitudes via 
Eqs. p3|) and Apparently, both contributions are 

very small. In addition and independent of the result 
that low-momentum 3N interactions are perturbative for 



cutoffs A < 2fm 



we find here that the contribu- 



tions of 3NFs decrease rapidly with increasing rank of 
the normal-ordered terms. 

The small contribution from the residual three-body 
Hamiltonian is the most important result of our study. 
It suggests that one can neglect the residual terms 
of the 3NF when computing binding energies of light 
nuclei. This is not unexpected and has been antici- 
pated in several earlier studies. Mihaila and Heisen- 
berg [III computed the charge form factor for ^^O within 
coupled-cluster theory and found a very good agreement 
with experimental data by considering only the density- 
dependent one- and two-body parts of 3NFs. Similarly, 
Navratil and Ormand [s^ observed in no-core shell- 
model calculations that density-dependent two-body 
terms are the most significant contributions of effective 
three-body forces. Our finding also support Zuker's [l^l 
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FIG. 8: (Color online) CCSD(T) results for the binding en- 
ergy of ''He as a function of the oscillator spacing and for 
model spaces consisting of A'' = 3 to A*' = 6 oscillator shells. 
The contributions from 3NFs are limited to the density- 
dependent zero-, one-, and two-body terms and exclude its 
residual three-body terms. 




N 



FIG. 9: (Color online) Data points: CCSD(T) results (taken 
at the huj minima) for the binding energy of *He with 3NFs 
as a function of the number of oscillator shells. Dashed lines: 
Exponential fit to the data and asymptote of the fit. Full line: 
Exact result. 



idea that monopole corrections to valence-shell interac- 
tions are due to the density-dependent terms of 3NFs. 
Note finally that the modeling of three-body interac- 
tions in terms of density-dependent two-body Hamiho- 
nians has a long history, see e.g. Ref. [6l[. Note that all 
these examples and the present study employ sufficiently 
"soft" or "effective" interactions. We expect that the 
smallness of residual 3NFs is a property of such interac- 
tions. We will study the cutoff dependence of this finding 
in future work. Finally, the smallness of residual 3NFs is 
also encouraging for future improved nuclear matter cal- 
culations, which currently include low-momentum 3NFs 
through density-dependent NN interactions [sJl. 

The smallness of the residual three-body terms is also 
for coupled-cluster calculations a most welcome result. 
This is attractive for two reasons. First, the inclusion of 
the residual three-nucleon Hamiltonian, as described in 
Subsection 111 Bl is computationally expensive. It exceeds 
the cost of a CCSD calculation for two-body Hamiltoni- 
ans by a factor of order 0(n„) -I- 0{nl) and is there- 
fore significant for a large number of unoccupied orbitals 
and/or large number of nucleons. Second, the omission 
of the residual three-body Hamiltonian will allow us to 
treat 3NFs within the standard coupled-cluster theory 
developed for two-body Hamiltonians (after normal or- 
dering). As a result, we can take the CCSD calculations 
one step further and include perturbative corrections of 
three-particle-three-hole clusters [62j . 

Let us neglect the residual 3NF terms of Eq. ^ and 
perform CCSD(T) calculations for the binding energy of 
*He. The approximate inclusion of three-particle-three- 
hole clusters improves the accuracy of our calculations. 
Our results are shown in Fig. [51 The comparison with 



the CCSD full lines in Fig. O shows that the triples cor- 
rections add about 100 — 200 keV of additional binding 
energy (at the minimum) for fixed number N of oscillator 
shells, and somewhat weakens the hio dependence. 

An exponential extrapolation of the (approximate) 
CCSD(T) minima to an infinite model space is shown 
in Fig. [5] and yields E^o = -28.24 MeV. This is in excel- 
lent agreement with the exact Faddeev-Yakubovsky re- 
sult E = —28.20(5) MeV. In our largest model space 
at the minimum Huj — 17 MeV, the ground-state ex- 
pectation values for the center-of-mass Hamiltonian is 
{Hem) ~ 20 keV while the expectaion value for the angu- 
lar momentum is zero for a closed-shell nucleus by con- 
struction. These results are very good, and it remains to 
be seen whether a more sophisticated treatment of triples 
excitations [i^l would lead to further improvements. We 
expect that the expectation value of the center-of-mass 
Hamiltonian decreases with increasing size of the model 
space. 

Finally, we also show the size of the CCSD(T) correc- 
tions in Fig.[7]as the horizontal dotted line. Clearly, these 
contributions are more important than the contributions 
from the residual 3NF terms, and this observation fully 
justifies the omission of the latter. 

IV. SUMMARY 

We have developed coupled-cluster theory for three- 
body Hamiltonians in the two-particle-two-hole cluster 
approximation (CCSD). We derived the corresponding 
coupled-cluster equations directly in a factorized form 
and thereby avoided the explicit construction and analy- 
sis of a considerable number of diagrams that enter these 
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equations. The resulting formulae were used for a very 
efficient numerical implementation. 

We have performed ab-initio coupled-cluster calcula- 
tions based on low-momentum NN and 3N interactions 
for the binding energy of ^Hc and compared to the ex- 
act Faddeev-Yakubovsky result. The 3NF contributions 
to the zero-, one-, and two-body terms of the normal- 
ordered Hamiltonian are dominant. The contributions 
from residual 3NFs are smaller than the corrections due 
to three-particle-three-hole cluster excitations and can 
therefore be safely neglected. Future work will include 
all 3N partial waves and studies of the cutoff dependence 
and of the convergence properties using low-momentum 
interactions with smooth cutoffs. Our findings tremen- 
dously simplify the computational cost of coupled-cluster 
theory with 3NFs. This opens the avenue to explore 
3NFs in medium-mass nuclei and to investigate questions 
related to modern nuclear interactions. 
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